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Abstract 

We propose a non-linear extension of the Fierz-Pauli mass for the graviton through a functional of the 
vielbein and an external Minkowski background. The functional generalizes the notion of the measure, 
since it reduces to a cosmological constant if the external background is formally sent to zero. Such 
a term and the explicit external background, emerge dynamically from a bi-gravity theory, having 
both a massless and a massive graviton in its spectrum, in a specific limit in which the massless 
mode decouples, while the massive one couples universally to matter. We investigate the massive 
theory using the Stiickelberg method and providing a ’t Hooft-Feynman gauge fixing in which the 
tensor, vector and scalar Stiickelberg fields decouple. We show that this model has the softest possible 
ultraviolet behavior which can be expected from any generic (Lorentz invariant) theory of massive 
gravity, namely that it becomes strong only at the scale A 3 = . 
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1 Introduction and Discussion 


Motivated by the observed accelerated expansion of the universe [1], and by the theoretical difficulties 
in ascribing it to a cosmological constant, there has been considerable activity in modifications of 
gravity at large scales in the past years. For instance, an accelerated expansion can be achieved 
in bi-gravity models [2], in models in which the Lorentz symmetry is broken by the gradient of a 
field [3], or in four dimensional models embedded in extra dimensions, as the self-accelerating DGP 
branch [4-6]. Some of these proposals have similar properties to massive gravity, which is probably 
the most straightforward and best studied modification of general relativity. 

At the linearized level, massive gravity is obtained by adding to the Einstein-Hilbert action a 
mass term for the metric perturbations The quadratic Lagrangian for this massive 

spin-two tensor field is given by [7] 

L = ^Ml,[h>^'^(h,^y + ahf,y-h^^,,ay-h^y,a^ + g^yhag,'^f^ -g^yah^+ml(h^-h^''h^y')] , (l) 

where h = g^'^h^y is the trace of the metric perturbation. There is a very stringent experimental 
bound on the graviton mass: nig < 7 x 10“^^ GeV [8], which is close to the inverse of the size of the 
observable universe. As already observed by Fierz and Pauli (FP) [7] the relative sign between the 
two mass terms is hxed uniquely by the requirement of having a ghost-free Lorentz-invariant (linear) 
theory. 0 The massless linear theory nig = 0 can be uniquely extended beyond quadratic order using 
the requirement of general covariance leading to the familiar Einstein-Hilbert action. But because the 
mass term breaks covariance, it has no unique non-linear extension. 

Govariance can be restored by introducing additional degrees of freedom, as for instance it is 
done with the Stiickelberg method [10]. Another approach is to introduce a second metric into the 
theory [2, 11, 12]. When one of the two metrics obtains a background expectation value, a mass 
term for the other metric is generated. Even though such bi-gravity theories are covariant, their 
completion of the Eierz-Pauli mass term is far from unique because one can write down an infinite set 
of invariant non-linear interactions between the two metrics. It is possible to obtain more uniquely 
defined bi-gravity theories: Ref. [12] considers a bi-gravity model described in terms of the vielbeins 
(tetrads), rather than metrics. Besides Einstein-Hilbert actions for both sectors it includes all possible 
cosmological constant like terms, that can be written down using these two vielbeins. The model [12], 
reviewed in appendix [Aj has two spin-two fields, one of which is massless, while the other has a mass 
term of the EP form ([T|). Interestingly, the model admits a limit in which the massless mode decouples, 
while the massive one couples universally to matter. 

The goal of the present paper is to investigate the resulting model of massive gravity. We would 
like to perform an analysis beyond the linearized level, and to compare our results with those obtained 
for generic massive gravity theories. In particular we want to investigate when does this massive 
gravity theory become strong? This question is important because it is related to the Van Dam- 
Veltman-Zhakarov (vDVZ) discontinuity [13]. The propagator of a massive graviton does not reduce 
to the massless one in the limit of vanishing graviton mass nig 0. The discontinuity between 
these propagators results in a discontinuity between the perturbative interactions of the massless 
and massive theories. However, precisely because the interactions of massive gravity become strong, 
this does not necessarily mean a discontinuity between the final nonperturbative results [10,14,15]. 
Ref. [10] showed that the scale at which any non-linear completion of the FP mass term ([T]) becomes 

richer structure of ghost-free mass terms is possible if one is willing to give up Lorentz invariance [9]. 
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strong never exceeds A 3 = {rin?gMpY^^ . Any generic completion, which becomes strong at a smaller 
energy, can be improved by adding suitable terms, to result in a theory which becomes strong at A 3 . 
This procedure is in general rather involved. Quite remarkably, we will prove that our model becomes 
strong precisely at the scale A 3 , without the need of such additional terms. 

The plan of this paper is the following. Section [5] describes the model of massive gravity we want 
to study in this paper using the vielbein formalism. In Section [3l the Stiickelberg analysis of [10] is 
extended to the case of the vielbein. In general, the vielbein formulation of gravity contains more (non- 
dynamical) fields than the standard metric formulation, which are compensated by additional gauge 
symmetries (the local Lorentz transformations). We provide the (unique) transformation which fixes 
these additional degrees of freedom. In Section U] we determine the quadratic action for the graviton 
and the Stiickelberg fields. We present a gauge choice, that explicitly decouples the scalar, vector and 
spin-two degrees of freedom. This choice, which, to our knowledge, had not been provided so far in 
the literature, leads to particularly simple propagators for the different polarizations. In Section [5] 
we classify the dominant interaction terms, and we show that the scale at which the theory becomes 
strong is A 3 . In that Section we also compute the tree level amplitude of the 2^2 scalar scattering, 
because it can be considered as a typical diagram used to compute the scale at which massive gravity 
becomes strong. However, we show, by a specific choice of parameters, that this interaction can be 
made to vanish at this scale, while other interactions still remain strong. These results are summarized 
in the concluding section [ 6 j In Appendix Owe explain how the model of massive gravity studied in 
this paper can be obtained via a decoupling of the bi-gravity model introduced in [12]. The following 
appendices are more technical. Appendix [B] describes some useful properties of a special product 
that generalizes the dehnition of the determinant, which is used to write the interaction between the 
two initial gravitational sectors. The Appendices [O and [D] contain useful intermediate results for the 
computation of the action in terms of the Stiickelberg fields. 

2 A massive gravity theory described in the vielbein formalism 

In this Section we outline the model of massive gravity that will be studied in this paper. The action 
of this model is 

S = j d‘^x^^/^R + 3ml(^{e - 7/)^(e + , (2) 

where we have denoted by the metric associated to the vielbein . The theory describes a 
massive graviton and is characterized by the Planck mass Mp and the graviton mass rUg. To write 
the cosmological constant-like interaction term we have introduced the notation 

(ABCD) = - A^aBpbC^cDsd , (3) 

with being the totally anti-symmetric Levi-Civita tensor. (Some useful properties of this product 
are collected in appendixjBj) This coupling term is a generalization of the measure, since {A^) = | A| = 
det(A). Using these properties it is not hard to show that the action can be rewritten as 

*5 = J {b + 3mfj - {[ef - [e2])| . (4) 
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where [e”'] is the trace of {rj ^ e)” . Hence, the hnal term can be interpreted as a FP mass term for 
the vielbein. By formally inverting the definition 

QfJ.u — (^un ) ( 5 ) 

we could have equivalently expressed it in terms of the metric. However, the resulting expression 
would appear as a complicated and unmotivated power series of the metric. This model is obtained by 
considering a specific limit of the bi-gravity theory introduced in [12] as is explained in appendix 
The specihc choice of the action (I2|) is obtained by imposing the reflection symmetry e^m — > , 

and by requiring that the Minkowski background is a solution. Indeed, from this last 

requirement, we see that the second term must factorize the (e — •rff' combination, and then the 
second factor simply follows from the reflection symmetry. This specihes the action ([2]) uniquely 
among the several interaction terms which may constructed starting from the vielbein, the Minkowski 
background, and the product (|3]). 

However, there is still a large arbitrariness in this procedure, given by the choice of the background 
term. In the action ([2]) we chose it to be equal to the Minkowski background in Cartesian 
coordinates. But, already choosing it to be the Minkowski metric in spherical coordinates would 
result in a different theory, since this term explicitly breaks covariance. We can also obtain more 
general solutions, starting from a different background metric in the action ([2]), and adding the 
corresponding source term which would lead to that solution in the standard case since the “mass term” 
is quadratic in it would not affect the validity of the solution. However, for definiteness, 

in the remainder of this paper we will only consider the case where the background is Minkowskian 
using Cartesian coordinates. 

The interaction term in dJ]) constitutes a particular completion of the FP mass term. Inserting the 
expression 


— 




( 6 ) 


in the second term of (l2|) gives 

AS = [iff - 1/2]) + i([/]3 - 3|/][/2] + 2 1/2]) 

+ ^([/l‘ - + 'Hf? + 8[/llfl - 6[/‘]) } . (7) 

We stress that this expression is exact (rather than just a perturbative expansion to fourth order). In 
the remainder of the paper we discuss the effects of the nonlinear interactions in setting the scale at 
which the model becomes strong. 

3 Stiickelberg fields 

We investigate the graviton interactions in the model outlined in the previous Section. Nonlinear in¬ 
teractions of massive gravity can be most easily studied through the Stiickelberg formalism, developed 
in ref. [10]. To do so for our model, we first have to formulate the Stiickelberg formalism in terms of 
the vielbein, rather than the metric (as done in [10]). In this Section we describe this computation in 
some detail. 
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The Stiickelberg formalism consists in performing a series of transformations in the action (j^ . 
and in promoting the parameters of these transformations to new fields. These fields appear in the 
new action together with additional symmetries, so that the original action can be recovered with a 
particular gauge choice. However, we can also choose alternative gauges in the new action, where the 
nonlinear interactions can be computed more easily. 

More specifically, we start from a symmetric vielbein perturbation and perform a general coor¬ 
dinate transformation x ^ y = y{x) combined with a local Lorentz transformation, with Lr]~^L^ = rj. 
Performing these transformations into the action ([2|) results in the replacements 


d^x 


d^y 




o ,, ^ab y ^ca ■ 

oyi^ 


( 8 ) 


We note that, even though e^m is assumed to be symmetric, is certainly not automatically 
symmetric. While the first term of ([2]) is invariant under these combined transformations, the whole 
action transforms into the Stiickelberg form 


S. 


st 



—gR + 



r]f{e' + , 


(9) 


where has been defined in ([S]). Notice that ([5D introduces the Jacobian in front of the cosmological 
constant term. 

The parameters in the local Lorentz transformations are not dynamical, therefore they can be 
integrated out using their algebraic equations of motion. Writing 


Lmn{x) = Lmpix) exp(7? ^b{x)) 


IP 


( 10 ) 


where is anti-symmetric, we expand the action ([9]) to first order in b: 


6bSst = GmlMl. 


j d'^x ^ (e'^ [e'rj ^6) — rfe' [e'rj ^6)^ . 


( 11 ) 


Using the property (IB.61) . it is immediate to verify that the first term in (1111) vanishes, due to the 
antisymmetry of bmn ■ Using (|B.8I) . the remaining term rewrites 

6bSst = - ^mlMp J d‘^x , ( 12 ) 


where we have used the notation [e'] for the trace, see below ([7|). This contribution vanishes if e' 
is symmetric, and therefore L = L is an on-shell solution. This symmetry of e' together with the 
requirement that L is a local Lorentz transformation, 

= ^er]~^L , and Lr]~^L^ = g, (13) 

determines L uniquely in terms of the other fields. See AppendixOfor a perturbative construction of 
L. The freedom in the choice of its sign is fixed by requiring, that L = 1 if {dx/dy)e is symmetric. 

Therefore, after we solve eqs. (fT3]) for L, and we substitute the solution back into Q, we are left 
with an action which explicitly depends only on the three dynamical fields a^, and (j). The latter 
two fields are obtained by decomposing 


y^{x) = x^ + a^{x) + d^(j){x) , 


(14) 
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into the spin-zero, (j ), and spin-one, , polarizations of the massive graviton (this decomposition 
introduces an additional 17(1) gauge symmetry in the theory). These fields are the starting point for 
the Stiickelberg analysis of the nonlinear interactions performed in the next Sections. 

To summarize the Stiickelberg formalism we have employed in this Section: In general, there are 
more degrees of freedom in the vielbein than in the metric, which are compensated by the local Lorentz 
transformations. This means that in the Stiickelberg description additional fields are introduced for 
both the general coordinate transformations and the local Lorentz transformation. But contrary to 
the Stiickelberg fields associated with the general coordinate transformations, the ones of the local 
Lorentz transformations are always auxiliary (i.e. non-dynamical) fields. When we enforce their field 
equations, we ensure that the composite vielbein e ', defined in ([8|), remains symmetric also after the 
Stiickelberg fields and 4> have been introduced. Hence the number of physical degrees of freedom 
in the vielbein formulation is the same as in the metric formulation discussed in [10]. 


4 Graviton and Stiickelberg propagators 


The aim of this Section is to compute the propagators of the massive gravity theory defined in Section [2] 
using the Stiickelberg field decomposition discussed in the previous Section. As usual the propagators 
can be read off from the quadratic action in the perturbations. However, we will encounter a few 
complications: First of all, the scalar cj) obtains a regular kinetic term only after a Weyl rescaling 
of the graviton field [10]. A second difficulty is that, before inverting the kinetic operator, we need 
to fix the general coordinate covariance and the additional 17(1) gauge symmetry generated by the 
Stiickelberg procedure. Even after the Weyl rescaling, the tensor the vector and the scalar cp mix 
with each other at the quadratic level (so that they cannot be used as such to describe independently 
propagating degrees of freedom). In the following we show how the last two problems can be solved 
together by choosing’t Hooft-Feynman-like gauge fixing terms. 

At quadratic order in f^i,, , (p the action Q is computed using (IB. 81) and the second order 

expansion of the matrix L~^ given in (jC.16|) . The result takes the rather complicated form 


52 =^Mf, I (fx 




r Uu + n/M. - r 




- r 




+ - °/) 


(15) 


+ 


f - + 2rd^d,(P - 2 /□</.] | 


The sign of the kinetic term of the vector is the standard one (had the sign of the mass term in ([2]) 
been opposite, would have been a ghost). The last two terms on the second line of (fT^ are the 
only two in which the scalar (p appears, and are not regular kinetic terms for a scalar. Following [10] 
we perform the linearized Weyl rescaling 

ffMu = Uu - i^rnlr]f,^(p (16) 
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of the graviton, to obtain a regular kinetic term for the scalar (p. The quadratic action (11511 becomes 


S 2 = -Mj. 


id‘^x 


r 




, “/t 


- □ 


+ 


f - r"// 




+ m„ 


{d^au + dya^) - 2 - 


+ rr?g\p + 


• (17) 


Even though we now have obtained a regular kinetic term for the scalar p, this action is still rather 
complicated and contains quadratic interactions between the three fields f^y , , and p. 

In addition, this action (1171) is invariant under the (linearized) general coordinate transformations 
and a t/(l) gauge symmetry 

Sffiu = ^ - ^rnlij^yil) , 5a^ = , Sp = - ip . (18) 

These gauge symmetries can be fixed by using the following gauge hxing functionals: 

GCC : = r,y-\f,^ , U{1) : © = . (19) 

Indeed, under the combined gauge transformations (I18p these functionals transform as 

50^ = i (□ — w?g) , (50 = (□ — m^) . ( 20 ) 

Note that at this order 0^ only transforms under general coordinate transformations, while 0 only 
under the 17(1) symmetry. This shows that by suitable gauge transformations these gauge hxing 
functionals can be made equal to any prescribed functions. When one performs an analysis at the 
level of equations of motions, it is most convenient to simply set both functionals 0 ^ and 0 to zero. 
However, since in this Section the aim is to obtain simple forms for the propagators, we employ the 
gauge hxing functionals (fTOll to dehne the gauge hxing action 

Sgf = - Mj,j d"x{0'^0^ + i®'} • 

The use of these gauge hxing functionals can be viewed as a generalization of the’t Hooft gauges in 
spontaneously broken gauge theories. In principle, we can allow for arbitrary normalizations in front 
of each of the terms in (1191) : As long as they are suitably chosen the tensor, vector and scalar modes 
remain decoupled at the quadratic level. However, in this ’’Feynman” gauge, the gauge hxed action 
reduces to the simplest form 



Hence in the gauge dehned by (f2T]) the different spin states f^y , and p are decoupled, and all have 
the same mass. In this gauge all components of these helds are dynamical. 

This action is very convenient: It is immediate from (I22p that the held redehnitions 


ffiu — fc^u ) 


mgMp 


3 m^gMp ’ 


(23) 
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3 point 

3 point 

scale 

{d‘^ (da) 

{d‘^py^{da) 

A4 

{d‘^py^{da) 

{d‘^p){da)‘^ 

An 1 
*^2 

{d‘^py{da) 

(dWf 

Ao 1 

2 

{d‘^py^{da) 

{daf 

As 

{d‘^py{da) 

{d‘^p){da)f 

^3 



As 


{d‘^p){daY 

As 

{d‘^p){da)‘^ 

{d‘^p){da)‘^ 

As 


3 point 3 point 

scale 

{d‘‘^py^ 

As 

{d'^pY {d‘^p)‘^{da) 

A41 

2 

{d'^p)^ {d‘^p){da)‘^ 

A4 

{d'^pf {d'^Pff 

A4 

{d'^p)^ (da)^ 

Ao 1 
'^2 

{d'^pf {d‘^p){da)f 

Aq 1 

2 

{d'^p)^ {da^f 

As 

{d'^pf {d‘^P)f 

As 


4 point 

scale 

{d'^PY 

{d'^pY{da) 

{d^YidaY 

{d’^pff 

A4 

An 1 
■^2 

As 

As 


[Ha mb 


Table 1: This table considers all possible interactions that give rise to four-point scatterings that becomes 
strong at scale A 3 or below, using canonically rescaled Stiickelberg fields (1^ . Such scatterings either arise 
directly from single four-point functions, ma^ or by combining two three-point functions, mb, with a scalar, 
vector or graviton propagator in between the two vertices. 


normalizes the fields canonically. Also the propagators for the graviton the vector field 

and the scalar A become very simple 

1 z 

^af)pa — TT {jlapTIpa T Vaallf^p TjalSVpa^ ^ j ^pu — Vpi '^ ~ 2 1 2 ' 

Z P ~r rUg 

These propagators can be obtained from massless propagators in the Feynman gauge by the straight¬ 
forward substitution: p^ + rri^. Even though the determination of the propagators here was 

performed in the vielbein formulation, it is obvious that it extends immediately to metric perturba¬ 
tions as well, since at the linearized level one simply has = 2fgy . Thus, this computation 
provides the gauge fixing, which was not explicitly given in ref. [ 10 ]. 

5 Dominant interactions at high energies 

After the discussion of the quadratic action and propagators, we now turn to the interactions of this 
massive gravity theory. The interactions of any massive gravity theory are rather involved, because 
there are many of them, and they all possess a rather complicated tensorial structures. Moreover, the 
interactions have a polynomial momentum dependence, and become non-perturbative at some “large” 
energy scale (in this context, “large” is defined with respect to the graviton mass). Therefore it is 
useful to classify which interactions are the dominant ones at high energies. 

As we reviewed in Section [3l in the Stiickelberg formalism the vector enters in the “pion” field 
with one derivative, while the scalar (p with two derivatives. Therefore, barring cancellations, the 
scalar couplings will be in general the largest at high energies. As was noted in [ 10 ], the (jxp (jxp 
scattering generally becomes strong at the scale A 5 = (rrigMp)^^^ . This scale falls into the regular 
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pattern of scales defined by 


Ap = (mP , (25) 

where p can be integer or half integer (notice that Ap is a decreasing function of p). All interactions 
that become strong at scale A 3 or below are grouped in table [H In particular, the scale A 5 appears 
in (pcj) —> (jicf) scattering by combining two three-scalar interactions. It was also shown in [10] that the 
scale at which this scattering becomes strong can at most be raised to A 3 . For a generic model of 
massive gravity, they sketched a procedure to obtain a suitable set of counter terms. We show that 
the model we are considering automatically satishes this property. 

To give a detailed discussion, we divide the presentation in three Subsections. The first one proves 
that our model does not have any interactions with either only scalars or at most one vector. In the 
second Subsection we show that all tree level processes become strong at a scale which is greater or 
equal to A 3 . In the final Subsection we explicitly compute (jicj) (jxf) scattering, as an example, and 
show that it is possible to extend the model such that this amplitude vanishes altogether. 

5.1 Absence of interactions with only scalars or at most one vector 

In the next two Subsections we want to show that our theory does not become strong below the scale 
A 3 . Here we prove that the most dangerous interactions, which have the schematic forms 
and (5^(/>)"■ (So), are all absent in our model. In these subsections we consider the scalar interactions 
before Weyl rescaling: The additional scalar interactions introduced by the Weyl rescaling will have 
an additional factor from (I16p . and have therefore the same strength as tensor interactions. The 
interaction term of ([91) can be written as 

^Sst = \mlMl I d^x(^{e - n)"(e + H)') , (26) 

where we have introduced the short-hand notation 

(27) 

where is expressed by ()14l) in terms of Op and (j). In the following, it is often useful to be able to 
resort to matrix notation to suppress the indices. Aside from the Minkowski metric 77, its perturbation 
/ and the Lorentz transformation L (which are matrices by definition), we encode the derivatives of 
y into the matrices 

where <I>p,y = and Fp,y = is the U{1) field strength. (For us vr is 

defined strictly by the matrix equation (1281) : our definition differs slightly from the conventions used 
in ref. [ 10 ].) 

Now let us first consider interactions with only scalar fields. To single out from (1261) the interactions 
((9^0)”, we can simply replace Hpj/ —> rj^u + and set both Lp,y and ep,y equal to the Minkowski 



metric (ignoring the vector and tensor contribntions). We can take L^i, = becanse in this 
case e' is antomatically symmetric. This gives 

ASst A J d'^x{4:r]‘^ + 477 $^ + $^) , (29) 

where we have employed the matrix notation defined in (j28p . All these terms vanish upon partial 
integration: For example, for the first term we obtain 

J d‘^x{4r]^<^^) = - J d'^xe‘^^'^eab^'' (t),cp4>,d<^ = j d'^x (I),p4>,cda = 0, (30) 

where we have used the dehnition (l3|) and the anti-symmetry of the tensor. Similar arguments 
also apply to the other two terms. Hence, all the interactions of the form (9^(/>)"■ vanish. Equivalently 
this result can be obtained by going to momentum space and realizing that then all matrices are of 
the form ()B.10p for which the angular bracket vanishes as is proven in Appendix [Bj As anticipated 
below (I7|), this shows that all these pure scalar higher derivative interactions naturally vanish, which 
was the motivation in ref. [18] to consider these combinations. 

Also the interactions of the schematic form (5^(/>)"■ (9a) vanish. To see this, we employ the matrix 
notation (|28p , and we again set e = 77 in the interaction term (|26p . We first observe that all the terms 
which are linear in F vanish. This is due to the fact, that all the other tensorial structures which 
would multiply F, namely 77 or 4>, are symmetric, while F itself is anti-symmetric. Therefore, we can 
set both and F^^ to zero, without losing the term that we are looking for. Doing so, we have a 
symmetric pion matrix, vr = . Once we insert it in eq. (I28p . we obtain a symmetric 

e' vielbein already, therefore L = 77 . The expression we are looking for can be found by evaluating 
eq. pHP for F^y = = 0 , L^j^y = rj^y . This results in an interaction action which is again of the form 

(1291) . upon the replacement <h —> ■ The final step is to single out from this expression 

the terms which are linear in the vector field. It is clear that they are of the form . Using a 

similar partial integration procedure as for the terms (<I>” 77 ^“") presented above (or again using (jB.lOp 
in momentum space), it is easy to verify that also such terms vanish. 


5.2 All interactions become strong at scale A 3 or above 

We now show that there are no scattering amplitudes that become strong at a smaller scale than A 3 . 
We only analyze the theory at the classical level, in particular, we do not consider loop graphs. We 
first consider the S-matrix elements that correspond to diagrams that contain a single vertex Vng,na,n^, 
with Ug tensor, Ua vector and n,f, scalar external legs. El Taking into account that at a vertex there is 
momentum conservation encoded in a single overall four dimensional momentum delta function, which 
scales as 1 / E'^ (where E is the energy in the scattering), we see that such vertex scales as 








/ E \ 

f E‘^ ’!•</. 

\Mp) 

\mgMp) 

\mjMp) 


(31) 


where the E in parenthesis arise from derivatives (scalars enters in the Stiickelberg formalism with 
two derivatives, while vectors with one), and the denominators arise from the canonical normalization. 

®We would like to remind the reader that, for reasons explained in the previous subsection, we count the scalar and 
tensor legs before Weyl rescaling. 
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To understand at which energy E such process becomes strong we need to compute the scaling of the 
corresponding S-matrix element Sng,na,n^, in which we integrate over all possible external momenta. 
Because the external particles are all on-shell, we integrate over 


J d^Vi^ipf -mf) = j ~ , (32) 

where in the last step we used that in the high energy regime we are considering E » rUg . Hence 
the corresponding S'-matrix element scales like 



(For notational convenience we have taken the square root of the 5-matrix element; the interaction 
becomes strong when the matrix element, or its square root, exceed one). 

Now we analyze the strength of general interactions with a single vertex. Because we are interested 
in scatterings, the total number of external legs Ug + Ua + > ?>. In total we can distinguish five 
different cases: (i) ria = Ug = 0 , (ii) Ua = , (iii) ng,n(j, > 1 , (iv) Ua = 1 and finally (v) 

na > ■ We have shown in the previous Subsection, that the model does not have vertices with only 

scalars, hence there are no processes that correspond to the first case (i). The second case (ii) involves 
interactions among only tensors. The corresponding 5-matrix elements 


1/2 /mgY / E \rig-4 

ng,na=n^=0 ~ \MpJ \MpJ 


(34) 


are bounded from above for < 4, and therefore do not lead to strong coupling behavior. For 
Hg > 5 these interactions become strong at an energy scale, which is greater than the Planck scale 
Mp. Next, let us consider interactions of type (iii) which involves at least one scalar and one graviton. 
By rewriting the general expression (j33p as 


1/2 / E \ng-l / E \2na / E \3(n,^-l) 

~ 

we conclude that these interactions are still weak at energy scales lower than A 3 , because all the ex¬ 
ponents of the factors are nonnegative, and the scales A 3 < A 2 < Ai = Mp are ordered hierarchically, 
see their dehnitions (1251) . Let us turn to interactions of type (iv), with a single vector, Ua = 1. From 
the previous Subsection we know that if there are no tensors, Ug = 0, the amplitude vanishes. The 
situation with at least one tensor and one scalar constitutes a special case of (iii). The remaining pos¬ 
sibility of type (iv), no scalars and rig > 2, becomes strong above the Planck scale Mp (analogously 
to case (ii)). This is clear when we rewrite 



1/2 /mg\2 / E \ng-2 

ng>2,na=l,n^=0 \MpJ \MpJ 

Finally, also case (v) remains weak at all scales lower than A 3 , because we can write: 


(36) 


0I/2 

‘^ng,na>2,n^ 


E \ng f E \2{na-2) / E \3n^ 

a7/ va^/ va^/ 


(37) 
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Figure 1: This diagram displays the three diagrams that contribute to the four scalar scattering at the scale 
A 3 , which result from (1551) . In the first diagram the scalar (j) is exchanged, while in the second diagram the 
graviton is the mediating particle. The last diagram results from the four point interactions of scalars. 


Hence, we showed that, in this model, all n—leg interactions with a single vertex become strong at 
the scale A 3 or higher. 

Finally, let us discuss tree diagrams that contain more than one vertex. Such diagrams can be 
obtained recursively by combining tree level diagrams with less vertices inside. Whenever we combine 
two such diagrams, we loose two external lines, and hence two factors of E'^ in the S'-matrix. At 
the same time we gain a factor E'^ in the amplitude, since we have an additional momentum integral 
over a propagator: f d^p/p^ ~ E"^ . (Because we consider tree diagrams all momenta inside a given 
diagram are fixed by momentum conservation.) Therefore, the scaling of the S'-matrix element of 
the combined diagram is the same as of the original two disconnected diagrams. We can repeat this 
argument recursively, as we split any tree level diagram in a series of single vertex diagrams. Since we 
already saw that all single vertex scatterings become strong at least at the scale A 3 , this argument 
shows that this is the case for any tree level scatterings. @ 

5.3 Can the (f)cf) — (f)cf) scattering amplitude vanish? 

The interaction vertices given in Subsection 15.11 can be employed to obtain several four-point ampli¬ 
tudes which (if not vanishing) all become strong at the scale A 3 . These amplitudes correspond to 
the scattering processes (/xp —>■ (fxj), (j)f —>■ (pf, acj) acj), aa —>■ aa, and cpf ^ aa, plus the crossed 
processes. As an example, we compute the (pcf) (fxf) scattering at tree level. In particular, we want to 
investigate whether it is possible to have a model where the full leading (pep (pep scattering vanishes 
at tree level. 

As will become clear below, the amplitude for epep (p<p does not vanish in the model ([ 2 |). Therefore, 
we consider the slightly generalized interaction term 

j d^x</(e-??)^(e2 + {2 +a)rje + {1 + , (38) 

where the real parameters a and f3 are arbitrary. The normalization factor 4 + a + P in (I38j) is chosen 
such that rrig still represents the graviton mass. Note that for a = P = 0 we recover the terms in 
eq. ([ 2 ]), which is symmetric under reflection of the vielbein. (Also in this more general model, the 
interactions which can potentially become strong at a scale lower than A 3 vanish for the same reasons 
we have discussed in the previous Subsections.) To be able to directly compare various scattering 

®A11 the arguments in this Section ignore any possible interference between different diagrams. Interferences can only 
soften the scattering amplitudes, so they do not affect the conclusion that the interactions are weak at energies smaller 
than A 3 . 
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amplitudes, we use the canonically normalized fields defined in (1231) . The leading interactions of the 
action (1381) can be expressed as 


AL 


gen 


D - 


—--— I 

m? Mp I 


2 + 36 . 2 ^ 1 1 + 36 


[FcFc] 


+ 


1 1 + 36 (J_ 

3 m?gMp 1 ^/6 


- |$S) - + (AUS - 1 (/J M - [/Ala<^=} 


+ 


a — 2b 
mjMj, 


+ 6D0,[Ad>2]) + i([ch2][F2] _ 


- 3(n0j2 ^ 3 [^2] Jj ^ 2 [/,] [ch3] - 6 [h^l] 

(□0,)2[F,]2) - _ 3n<^^[cl,2] + 2[ch3])} 


where we have defined the parameters 


a 

4 + a + /? ’ 


6 
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4 + a + /3 ' 


(40) 


To obtain these interactions we first performed a Stiickelberg transformation to the action (|38p . we 
then expanded it up to first order in / and fourth order in vr, and we finally performed the linearized 
Weyl rescaling. To single out only the leading terms, we then substituted tt = A + ^ and kept only 
those terms with the highest power of 4>, since they contain the greatest number of derivatives. In 
this way, we found the expansion (|C.18p given in Appendix [C] for IT defined in (I27p . Finally we worked 
out the brackets (...) in terms of traces [...], using the identities of Appendix iBl 

This fourth order leading expression is rather involved for generic values of a and 6 . However, 
notice that if a = 26 the entire third and fourth lines of (j39j) vanish. In particular, there are no pure 
scalar interactions anymore. This is also the case for the model ([2]), with o = 6 = 0, which we are 
mostly concerned with in this paper. However, if in addition 6 = —1/3, also the second line and the 
last term on the last line vanishes, and only three different interactions survive: one two-scalar vector 
interaction and two two-scalar two-vector interactions. In particular there are no interactions with 
tensors left. Hence the only possible four-point scatterings involve two scalars and two vectors. 

We compute the leading tree level tpcf) (pcf) scattering which becomes strong at the scale A 3 
for generic values of a and 6 , using the leading expansion (f39P . The process is described by the 
three diagrams shown in figure [TJ The computation of these diagrams is a straightforward exercise in 
Feynman graph computations and therefore only the result is given here. We neglect the graviton mass 
against the external and internal momenta of the scattering. In terms of the standard Mandelstam 
variables s, t and u, the amplitude reads 


M{cj)cj) (/(/) 


(1 + 36)^ ^ (1 + 36)^ , a-26 
8 12 ^ 3~ 


stu 

mjMj, 


(41) 


The three terms correspond to the diagrams [TJa, [TJb, and[TJc, given in figuredl respectively. The total 
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amplitude expressed in terms of the center of mass energy E equals 


A4{(j)(j) 4>(j)) = — - 8(2 b — a) + {1 + 3b) 




- M^nig 


sm 


9 


(42) 


where 9 dehnes the angle between the momenta of an ingoing and an outgoing particle in the center 
of mass frame. We see that, besides the special point a = 2b = —2/3, there is a whole parabola 
a{b) = 26 + (1 + 36)^/8 such that the full tree level (jxp (pcj) scattering vanishes at the scale A 3 . 
We emphasize that this four-point scalar scattering alone is not sufficient to understand the strong 
coupling dynamics of massive gravity theories, and of our model in particular. 


6 Conclusions 

The model of massive gravity considered in this paper is standard general relativity with a cosmological 
constant and a Fierz-Pauli mass term for the vielbein. By requiring a reflection symmetry of the 
vielbein, e^a —> this theory, like standard Fierz-Pauli massive gravity, is described by only 

two parameters: the Planck scale and the graviton mass. As such it constitutes one of the simplest 
non-linear extension of massive gravity theories. This model can be obtained from a bi-gravity 
theory with an extremely simple interaction term between the two sectors. The bi-gravity theory 
contains a massless and a massive graviton state, and it admits a limit in which the massless graviton 
decouples [12], as we show in appendix^ In the present work we studied the model of a single massive 
graviton that emerges in this limit. 

We have investigated whether the simplicity of the model has some physical implications, or only 
has aesthetic merit. To do so, we studied this model at the non-linear level (since, at the linear 
level, it is equal to the standard Fierz-Pauli model). As the present model is a special type of massive 
gravity, it shares the problems that massive gravity theories have in general. The main one is that they 
are plagued by ghosts at the the non-linear level [16,17], which can never be avoided at the quartic 
order in the perturbations [18]. A related difficulty is that the self-interactions of a massive graviton 
become strong at macroscopic distances from a source [14]. In general, the scattering becomes strong 
at the energy scale A 5 = (rrigMpY^^ ■ However, by adding to the starting theory a set of suitable 
nonlinear interactions, this scale can be raised up to A 3 = {m^Mp)^^^ at most [10]. Even though our 
theory does not solve these essential problems of massive gravity theory, it provides a minimal model 
that has self-interactions which become strong at the highest possible scale A 3 . Therefore, the model 
here discussed may be considered as a prototype of massive gravity, since, in addition to its minimal 
formulation, has the best behavior that we can hope to obtain for such theories. 

We conclude with a note on the use of the vielbein rather than the metric formulation. In eq. dZD, 
we wrote the (exact) lagrangian of the model in terms of metric perturbations. Ref. [18] already 
showed that special combinations of terms have the consequence of removing from the theory the self¬ 
interactions which involve only the scalar polarization of the graviton, before Weyl rescaling (such 
interactions - if present - would lower the strong scale below A 3 ). However, what in the metric 
computation appears only as a computational result (obtained by allowing for arbitrary coefficients, 
and then finding which combination eliminates the unwanted terms), in the vielbein formulation 
corresponds to one of the simplest combination of a Fierz-Pauli mass for vielbein perturbations and a 
cosmological constant. This suggests that, despite it is very rarely considered, the vielbein approach 
may be more suitable for the study of massive gravity, and, hopefully, for finding the improvements 
that it still requires. 
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Appendixes 

A Derivation from a bi—gravity theory 

In this appendix we explain how the massive gravity model discussed in this paper can be derived 
from the bi-gravity theory discussed in [ 12 ]. One starts from the two metrics 5 +^^ and which 

we rewrite in terms of two vielbeins and e_^m) using the standard definition 

g± limV ■ (-^d) 


The action for the model is 

Sbi = I {^Ml R+ + A+) + R_ + A_) - 2 Ao e^.)} . (A.2) 

The two gravitational sectors are characterized by the two ‘’Planck masses” M± and the cosmological 
constants Aq and A±. 0 The double covariance of the model is broken by the last term. 

Let us now proceed to the study of the spectrum of the model for the Minkowski background. We 
set 


6 ± fMU 


Tlfiy + f± 


fiiy ? 


(A.3) 


and we expand the action (|A.2p at the quadratic level in the perturbations f±^u- The resulting action 
is not diagonal in terms of these two modes; however, it can be diagonalized through the redefinition 


f+A = A 
f-,J r + iyi 



M. 


r = 






(A.4) 


The mode /o^i/ is massless, while has a Fierz-Pauli mass nig, which is related to the cosmological 
constant Aq by| 


Aq 



Mi 


Ml 


M+M_ 



(A.5) 


The Planck mass Mp (as obtained from the kinetic terms) is found to be identical for both gravitons. 
The graviton mass nig vanishes for Aq = 0, as a consequence of the enlarged covariance. 

The massless graviton decouples in the limit of either M_|_ ^ 00 (with finite M_), or M_ —> 00 
(with finite M_|_). Both the Planck and graviton mass appearing in (IA.5I) are finite in these limits, so 
that the quadratic actions for the massless and massive modes remain finite. However, all the nonlinear 
interactions involving the massless mode fogi, vanish in either limit. For instance, in the first limit, 
we have at leading order /_|_ oc /o/M_|_,/, and /_ oc /o/M+,//M^ (suppressing for shortness the 

we see below Ao needs to be positive to avoid a tachyonic mass for the graviton. 

®We correct a typo in the definition of the mass appearing in [12]. 
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tensorial indices). Therefore, writing the original action (IA.2I) in terms of / and /o, we see that any 
term containing the massless mode is suppressed by a negative power of M+, and therefore vanishes 
as M_|_ ^ oo. The only exception is the quadratic kinetic term from the + sector, which results in 
a finite kinetic term for the massless graviton with the Planck mass (lA.Sp . In this limit, the massive 
graviton is identified with /_, and a theory of a single massive graviton is obtained by coupling all 
matter fields in the — sector. Obviously, an analogous situation is obtained in the other limit. 

To obtain the model of this paper, we study the theory in either of the two (equivalent) limits, by 
ignoring the decoupled massless graviton. We also restrict our attention to the Minkowski background. 
For instance, in the M_|_ ^ oo limit, this gives , e-^jy = c~^e^y = c“^(? 7 ^iy + f^y), and 

the action (IA.2P reduces to ([2]). 

B Properties of (ABCD) 

In this Appendix we collect various helpful properties of the angular bracket {ABCD) defined in 
eq. m- First of all the ordering of the matrices A,... ,D is irrelevant in this expression, because of 
the two Levi-Civita tensors in its definition. For the same reason, this product is invariant under the 
simultaneous transposition of all four matrices. As we remarked in the main text, it generalizes the 
notion of a determinant, in the sense that (A'^) = det(A) = |A| . However, cannot be written 

as a determinant. The property of a determinant, that the determinant of a product of matrices is 
equal to the product of their determinants, generalizes to 

( {aAb) {aBb) {aCb) {aDb) ) = |a| {ABCD) |6| , (B.6) 

for any matrices a,b,A,B,C,D. 

The angular bracket can be rewritten in terms of traces [A] = rj^'^Aba ; the resulting expression is 
rather involved 

{ABCD) = i ([A] |B] [C] ID] - [A] |B] [CD] - [A] [D| [BC\ - \B] ]D] [AC] - [AB| [C] [D| 

- [A] [BD] [C\ - [AD] [B| [C] + [AB] ]CD] + \AC] [BD] + [AD| [BC] + [A][BCO] 

+ [A] [BDC] + [B] [ADC] + [B] [ACD] + [C] [ABD] + [C] [ADB] + [D] [ABC] 

+ [D] [ACB] - [ABCD] - [ADBC] - [ACDB] - [ABDC] - [ACBD] - [ADCB]j . (B.7) 

However, when some of its entries are equal to the Minkowski metric 77 , its expression simplifies 
considerably: 

(ry^) = 1, (ry^A) = i [A] , {rj^ AB) = 1 ([A] [H] - [AH]) , 

(ryAHC) = ^ ([A] [B] [C] - [A] [BC] - [B] [AC] - [C] [AH] + [AHC] + [ACH]) . (B.8) 

Also, we can use (IB.7h to express the determinant of a matrix in terms of traces: 

\A] = {A*) = P([A]- - 6]Af]A^] + 3]A^f + 8[A|[A=>] - 6[a"]) . (B.9) 
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Finally, if the matrices A,... D are formed from five arbitrary vectors p, q, r, s and t, 

Aaa — PaQa j ^fib — P/3^b j C^c — P-y^c ) ^Sd — PS^d j (B.IO) 

then, by the anti-symmetry with respect to the exchange of any two of the Pa,Pi3,P'y,P6 inside the 
bracket expression, we find {ABCD) = {rjABC) = {rf'AB) = 0. 

C Perturbative expansions 

This Appendix is devoted to some technical details of the perturbative expansions, that we use in the 
main part of the text to determine the interactions of the massive gravity theory in the Stiickelberg 
formulation. The interactions are encoded in the expression (1261) . where the matrix 11 is defined in 
eq. ()27l) . The first step for the computation of 11 is to determine L from eqs. (I13j) . Since the interactions 
of the graviton are determined by an expansion around the Minkowski background, we consider the 
infinitesimal general coordinate transformation 

F)'T' 

1 + (C.H) 

dy 

and we expand L = in a power series in e and its transpose (more accurately, is a sum of 

monomials of degree re, where each of the monomials is a product of e and its transpose). Using this 
expansion, eqs. (I13p can be rewritten as as two recursion relations 

n 

Ln+l — ~ , L^+i + ^ , (C.12) 

k=l 

which, altogether, determine L order by order in the expansion. From this condition, and from taking 
L = r] when the change of coordinate is trivial (e = 0), we Hnd the recursive solution 


1 * 

Lo = r] , Ln+i = - ep-^Ln - , (C.13) 

k=l 

where the last term in parenthesis must be evaluated only for re > 1. This determines L uniquely; up 
to cubic order, the explicit solution reads 

L = , - 1 (e - + 1 {3d - ed - d, - (df) 

+ i (e^e^ + cdt + ddf + dd - d-cd + (dfe + (df - bd) + ... . (C.14) 

In this expression the presence of between any two consecutive e or is understood. We can now 
use this information to determine the matrix IT = ry + Hi +112 + FIs +..., see (I27h , where the subscript 
indicates the order in which the graviton / (before Weyl rescaling) and the pion field tt appear in this 
expression. In fact, in this work we only need the expansion up to third order: We must evaluate the 
interaction term (|26p up to quartic order in the fields, since we are at most interested in four point 
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interaction vertices. However, since one of the two factors entering in (I26p does not have a background 
part, it is sufficient to expand H at cubic order in the fluctuations. 

From the definition (j27|) . we see that we need to invert L up to cubic order. The inversion up to 
cubic order in e is obtained in a straightforward manner from (|C.14p . However, one has then to realize 
that e itself is an expansion series in terms of the fields which are contained in / and vr. Therefore, 
we must now expand e at cubic order in the physical fields. This is done expanding the definition of 
e, eq. (1C.lip , where (in matrix notation) e = rj + f, while vr enters in the inverse of dx/dy as written 
in ([28p . The expansion up to cubic order for e reads 

e = / — vr — vrT/“^(/ — vr) + TTr]~^Trr]~^{f — vr) + ... . (C.15) 

and, finally, L~^ is found to be 

L~^ = r]~^ + ^ - vr^ + ^ (^3vr^ - (vr^)^ - vrvr^ - vr'^vr^ + ^ (^/vr^ + /vr - vr/ - vr^/^ 

+ — ^(vr"^)^ — 5vr^ + vr(vr^)^ + vr^vr^ + vrvr^vr — + vr^vr^ + (vr^)^vr^ 

2 vr^/ — /(vr^)^ — vr/vr + vr'^/vr^ — vr/vr^ — /vr^vr + vr'^/vr — /vr^ + /vrvr'^^ 

+ ^ + vr^/^) + ... . (C.16) 

In both this and the next expressions we have suppressed writing the matrix ri~^ between consecutive 
factors. Inserting this expression in eq. (1271) . we obtain 

n = ^ ^ ^ (^3vrvr^ - vr^ - vr^vr - (vr"^)^) + ^ (^/vr^ + /vr - vr/ - vr'^/^ 

+ — fvr^ + (vr^)^ — vr(vr^)^ — vr^vr"^ — vrvr'^vr — vr^vrvr^ + vr^vr^ + (vr'^)^vr') 

16 V / 

vr — /(vr^)^ + vr^/vr"^ + vr/vr^ — /vr^vr + vr^/vr — /vr^ + /vr-vr"^ — 2vrvr^/^ 

+ ^ (vr/^ - /'vr^ - /'vr + vr^/^) + ... . (C.17) 

which we finally insert in the interaction term (j26p . This gives the interactions between the various 
polarizations of the graviton for the model we are discussing. 

Eq. (13911 of the main text included the dominant interaction terms up to quartic order. These 
leading terms are obtained from the general expressions just given, by substituting (j28j) in the expres¬ 
sion above, and only keeping terms with at most a single / or two H’s. (All other terms have fewer 
derivatives, and so do not control the high energy limit of the model.) This gives 

H = + $ + i$(E -2 /) - ^(E -2 /)$ + ^(sAA^ - A^ - (A'^)^ - A^a) 

+ p p _ 2/$2 _ 2$v) • (C.18) 
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Notice that many possible structures are absent. In particular all terms with higher powers of and 
no other fields have canceled. The reason for this is that in (jC.171) terms that could give such terms 
have coefficients that add up to zero. For example, when substituting vr —> $ in the combination 
— tt^ — (vr^)^ — gives zero because <I> is symmetric. This completes the development of the 
expansions of the various functions that appear in the main text to the order required there. 


D Exact expressions 


In addition to the perturbative expansions presented in the previous Appendix, it is also possible to 
derive closed exact expressions. Such results can be obtained as follows: By multiplying the relations 


= e'L-y and . 


dy 


\dyJ 


we can obtain an equation for rj of which we can take the formal square root 


e = rj 


_i dx _i /dx\'^i 1/2 


Substituting this back into one of the equations (ID.191) , we find 

fdx\ 

J 




/dx\T'\ _^dx _i /dx\Ti-i/‘i 
V e( —) 


This in turn results in 


n = 


dy 


dx idy 


dx 


_i fdx\T _il-l/2(9x 
pA - 77 - 


■ ey e 


\dy) 


dy 


e . 


(D.19) 


(D.20) 


(D.21) 


(D.22) 
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